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1 Introduction 

One of the main aims of high energy nuclear colhsions is to explore a new state of matter, 
the Quark-Gluon Plasma(QGP), through heavy ion collisions in the laboratory. It was 
theoretically proposed that a suppression of j/ijj production in relativistic heavy ion 
collisions can serve as a clear signature for the formation of QGP [1]. Subsequently this 
suppression effect was observed by NA38 collaboration [2]. However, successive research 
pointed out that such suppression could also exist in hadronic matter (HM), even though 
caused by a completely different mechanism [3]. The anomalous J/^ suppression has 
recently been reported by the NA50 collaboration [4] and there have been a number 
of attempts to explain it [5-7]. Some authors believe that the data may implicate the 
possibility of the formation of a QGP [8]. For understanding the experimental data 
clearly, the dissociation mechanism of J/^ in hot QGP must also be studied carefully. 

In QGP, quarks and gluons are deconfined and the confining force between quark and 
antiquark vanishes, the only interaction between quark and antiquark is the Coulomb- 
type color interaction. The color charge of a quark will be screened by the quark sea in 
the plasma. Due to Debye screening the final yields of J /ip will be suppressed. Binding 
and dissociation of J/^ at a finite temperature have been studied in the non-relativistic 
formalisms [9,10]. The J/^ was regarded as a non-relativistic bound state in those papers. 
However, generally speaking, the motion of a quark and an antiquark in a meson is 
relativistic, even for charmonium. As pointed out in Ref.[ll], for the cc system, the 
kinetic energy is about 13% of the total energy and the ratio of relativistic corrections 
to the quark mass will not decrease with the increase of quark mass if the interaction 
is of Coulomb-type. As a result the bound state equation for J /ip in general should 
be relativistic. So it is an interesting task to discuss the binding and dissociation of 
charmonium in hot matter in a relativistic formalism. This is the main purpose of this 
paper. 

It is well known that the Bethe-Salpeter (BS) equation [12] is the only effective rela- 
tivistic equation of the two-body bound state problems. Because of its consistency with 
quantum field theory, the BS equation can be used for the study of the binding and dis- 
sociation of charmonium and bottonium which has seldom investigated in the framework 
of relativistic formalism in a thermal environment and of great interest. We expect that 
our calculations of observables could readily yield results at variance with those of non- 



relativistic models. Furthermore, we are not only interested in the energy spectrum of 
mesons, which is an important source to study the interquark dynamics, but also in the 
wave functions which play a key role in the calculations of the root-mean-square(r.m.s) 
radii of cc and bh bound states. 

In this paper, we shall discuss the binding and dissociation of heavy quark resonances 
in the hot matter within the context of the BS equation. In section 2, we focus our 
attention on the interaction between quark and antiquark in mesons and the properties 
of the BS equation. In section 3, we use the BS equation to calculate the mass spectra, 
r.m.s. radii and critical values of the screening masses for the cc and bb bound states and 
compare them with the previous results. The sensitivity of our results to the Lorentzian 
structure of the confining potential is also checked. In section 4, we discuss the results 
and conclude. 

2 Formalism 

It is well known that the BS equation is a proper tool for describing the relativistic 
two-body bound state problems[13]. The full bound state BS equation in momentum 
space, written in the two-sided notation, reads 

{riip + i) - mi)xp{p){ri2f - i> + T^2) = t:^ j d'^p'V{jp,p;P)xp{p') (1) 

where r/j = „^"1^„^. (i=l)2), Xp is the momentum-space wave function for the quark- 
antiquark system with total four momentum P in momentum space, p is the relative 
four momentum. V is the interaction kernel that acts on XP ^"^^ formal products of 
Vxp(p') in Eq.(l) take the form Vxp(p') = VsXp{p') +7^ ^l^yvXPijp')-, ^^ which Vg and 
Vy are scalar, vector potential respectively. The short-distance behaviour of the V can 
be calculated in QCD using perturbative theory. However, the long-distance behaviour 
of the V involves non-perturbative effects, and the Lorentzian structure of the confining 
potential is not established theoretically in QCD. Consequently, we shall treat the form 
of V in a partially phenomenological way. The parameters rrii, m2 should be interpreted 
as effective constituent masses and similarly the whole propagator is an effective one. 

Using the standard reduction and spin-independent treatment , one can get the spin- 
independent reduced Salpeter equation [14] for the three-dimensional equal-time BS wave 



function 

(^{p) = I dp\p{p^,p), (2) 



d^p' 



{M-E^- E^^ip) = j ^Y. Ft{p,p')V,{\p- p'UiA (3) 



(2vr)= 
where M is the mass of a qq bound state, Ei = {p^ + m'f)^ , i=l,2 represent a quark and 



l=S,V 



an antiquark, respectively. The functions F^* and F/* appearing in Eq.(3) are 



{p-p'f' Ei+mi E2+m2 Ei+mi -E2+^2 ^^ -. 

{E[+mi){E'^+m2) ^ E'^+m2 E[+mi E[+mi E'^+m2'^ ' ^ ^' 

(4) 

and 
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Here E[ = (p' + m'f)^ , i=l,2. Since -F^* and F?* are spin-independent, the singlet and 
triplet are degenerate. There is no coupling between different orbital angular momenta 
in Eq.(3). We can therefore extract the angular dependence of (j){p) in a single spherical 
harmonic basis 

Hp) = <PnL{\p\)YLm{h- (6) 

By using the following identity 

4tt 



2L + 



Y E ylm{P)yLn.{p') = PUcOSe), (7) 



where 9 is the angle between the p and p', one can obtain the following equation for the 
radial wave function (/>nL (|p|) 

{M^L - Fi - E2)cPnL{\p\) = J 1^ E Fr{p,p'm{\p-p'\)PL{cOsd)^,,U\p'\)- («) 



t=s,v 



Eq.(8) gives a well-defined eigenvalue problem for the masses of the qq bound states in 
momentum space. Here we would like to point out that the momentum dependence of 
the interaction is treated exactly in the above equation. 



In the non-relativistic limit, Eq.(8) can be reduced to the usual Schrodinger equation 
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{MnL - El - E2)cPnU\p\) = / J^^iVsiW - p'\) + Ki\p - p')]PUcOS0)<PnL{\p'\). (9) 



It will be convenient in calculating the r.ni.s radii of qq bound state to transform 
(prihilpl) to the position space. Making use of the following identity 

oo 

exp{ip-f) = iTTj2i'ji{pr)YC^,{e,^)Yi^,{e,ip), (10) 

1=0 

where the direction of jf and r is specified by the polar angles (G, ^) and (6, (p), respectively, 
we get 

4>nLm{r) = j^ J exp{ip- r)4)nLrn{p)d?'P 
= i^I P^dp JL{pr)<PnL{\p\)YLm{9, 4>)- 

where we have used the orthogonal relation of spherical harmonic function and Eq.(7). 
According to Eq.(lO) one can get the radial wave function in the position space 

-^"^(1^1) = ^2 / P^dpJL{pr)<PnL{\p\) (12) 

in which the special functions jo{x) = ^^, ji{x) = ^^ — ^^. Solving Eq.(8) and 
Eq.(12), one can obtain the masses and corresponding r.m.s radii of the bound states. 

To solve Eq.(8), one must have a good knowledge of the potential between quark 
and antiquark. At present, it is commonly accepted that the interaction between quark 
and antiquark consistes of a short-range part describing the one-gluon-exchange(OGE) 
potential and an infinitely rising long-range part responsible for the confinement of the 
quarks. As it is well known that the OGE potential is a pure vector interaction. However 
the Lorentzian structure of the confining interaction is not clear. Wilson loop technique 
suggests that the confining potential should be taken purely scalar [15], but relativistic 
potential calculations [16-17] showed a need for some vector confinement. Therefore, we 
choose a confining potential to be a mixture of a scalar and a vector [18]. This leads to 
the following potential 

V{r) = Vs{r) + K,(r), (13) 

with 

Vs = (1- x)a r, 

Vy = xa r - - — . 
6 r 



where a is the string tension, Oe// = — fos the effective coupHng constant and x the 
vector-scalar mixing parameter obey the condition < a; < 1. Note that the edge of the 
interval x = corresponds to the case of pure scalar confinement. 

In a thermodynamical environment of interacting light quarks and gluons at temper- 
ature T, quark binding becomes modified by color screening [1] 
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Here ^ is the Debye screening mass (which is assumed to be a function of temperature T) 
and the Debye screening length td is defined as the inverse of the screening mass, r£i = -. 
It is necessary to note that the factor exp(— ^u r) not only reflects the color screening effects 
but also avoids the infrared divergence. In fact, as pointed out in Ref.[19], the color 
screening effects are also required at T=0 in order to flt the experimental properties of 
quark systems. In this paper, we use the color screening potential (Eq.(14)) to study the 
binding and deconflnement of heavy quark resonances . In momentum space the potential 
can be written as 



Vs{\p-p'\) = {l-x) 
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(15) 



(16) 



The constants a, n, x and ag are the parameters characterizing the potential. 

3 Calculations and Results 

Based on the formula above, we first calculate the mass spectra and r.m.s radii of cc 
and bb bound states with the vector-scalar mixing parameter x = 0, which corresponds to 
the pure scalar confinement. The numerical results are listed in Table I. The data used 
in our studies consisted of the spin-averaged masses of bb and cc states and are given in 
the third column in Table I, with 

1 



Mnl 



4{2l + 



-^(2i + l)M(n,i,/, 



s . 



(17) 



Here we restrict ourselves to the first two radial excitations, corresponding to the (spin- 
averaged) J/V' and T for n=l, /=0, to the iJj' and T for n=2, /=0, and to the Xc and Xb for 
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n=2, 1=1. Therefore M„; can be written explicitly as: M„o = \[3M{n^Si) + M{n^So)], 
Mni = ^[5M(n3P2) + 3M(n3Pi) + ^M{n^Pi) + M{n^PQ)] (n=0,l corresponding to the 
first two radial excitations), where we recall the usual spectroscopic notation ^^^^Lj for 
a state with orbital angular momentum L, spin s, and total angular momentum J; S, 
P, ... correspond to orbital angular momentum L=0,1,..., respectively. Since the spin- 
singlets bb{^So) and bb(^Pi) have not yet been unambiguously confirmed by experiment, 
we have therefore used the results of previous spin-dependent fits to the data to estimate 
the centers of gravity of the incomplete multiplets [14]. However, Ref.[10] compared their 
numerical results with the spin-triplet rather than the spin-averaged masses of bb and cc 
systems. 



rable I. T 


le mass spect 


ra and rms radii of the cc and bb bound states. 






nl 


Data 


MruiGeV) 


< r^ >^(fm) 








Ref.[9] 


Ref.[10] 


Ours 


Ref.[9] 


Ref.[10] 


Ours 


cc 


IS 


J/V'(3.068) 


3.0697 


3.0700 


3.067 


0.4490 


0.4453 


0.2868 




2S 


V'' (3.663) 


3.6978 


3.6863 


3.663 


0.8655 


0.9034 


0.6317 




3S 


V'" (4.025) 


4.1696 


4.0806 


4.019 


1.2025 


1.3765 


0.8290 




IP 


Xc(3.525) 


3.5003 


3.5054 


3.526 


0.6890 


0.7000 


0.5144 


bb 


IS 


T(9.436) 


9.4450 


9.4310 


9.436 


0.2249 


0.2211 


0.1873 




2S 


T'(10.013) 


10.0040 


10.0083 


10.013 


0.5040 


0.4998 


0.4480 




3S 


T"(10.341) 


10.3547 


10.3564 


10.343 


0.7336 


0.7457 


0.6749 




IP 


Xb{9.899) 


9.8974 


9.8981 


9.901 


0.4041 


0.3982 


0.3569 



In our calculations, we have used the following parameters. a=0.22 GeV^, ^o=0.06 
GeV, mc=1.474 GeV, m6=4.762 GeV, as{cc)=0A7, and as{bb)=0.38. Ah these param- 
eters are within the scope of customary usage. According to the concept of a running 
gauge coupling constant, in the computations we allow different values of a^ for charmo- 
nium {as{cc)) and bottonium {as{bb)) [20]. The relative magnitude of as is in accordance 
with the idea of asymptotic freedom as expected for the strong gauge coupling constant 
of quantum chromodynamics, that is, as{bb)<as{cc). 

In order to compare our results with those given in Refs.[9-10], we calculate the quan- 
tity x^ which is defined as 



X 



2 _ l^nlV^J^nl ~ ^^^n 



nl 



(18) 



iV-1 
with N being the total number of {nl) state. We obtain x= 0.0066 GeV for our numerical 
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results. Comparing with x= 0.0223 GeV (for Ref.[10]) and x=0.0511 (for the Ref.[9]), 
one can observe that the mass spectra obtained at present are more consistent with the 
experiment than the previous results. 

As mentioned above, the wave functions play a key role in the calculation of r.m.s radii 
of cc and hh systems . The last column in Table I shows that the r.m.s radii of cc and hb 
bound states given by our calculations are smaller than those of Refs.[9,10]. This means 
that J/V' is more tightly bound in our case than estimated by non-relativistic models. 

Next, we study the dissociation of hb and cc systems. A suitable quantity to observe 
the vanishing of bound states is the dissociation energy 



a 



E^Ufi) = mi + ms + - - Mniip.). 



(19) 



The dissociation energy is positive for bound states and turns negative for the continuum. 
Thus 

ElLii^c) = (20) 

defines the critical value of fi beyond which there is no bound state for the given quantum 
numbers. The calculated results of the E'2lg{nc) for the cc and bb systems are given in Fig.l 
and Fig. 2 respectively. The figures show that the dissociation energies of our calculations 
(solid lines) are shifted to larger ^ regions in comparision with those (dotted lines) of 
Ref. [10]. The reason is that our calculations are based on the relativistic formula in which 
the momentum dependence is treated exactly, which is different from that used in Ref. [10]. 

The calculated critical values of Debye masses and Mni{pc) for cc and bb resonances 
are given in Table II. 

Table II. The calculated Hc and M^i for charmonium (cc) and bottonium {bb). 





state 


/Uc(GeV) 


M„K^c)(GeV) 






Ref. [9] 


Ref. [10] 


Ours 


Ref. [9] 


Ref. [10] 


Ours 


Charmonium 


IS 


0.700 


0.600 


0.900 


2.9145 


2.8779 


3.1911 




2S 


0.360 


0.260 


0.470 


3.1725 


3.2964 


3.4160 




IP 


0.342 


0.242 


0.450 


3.1982 


3.3513 


3.4363 


Bottonium 


IS 


1.560 


1.500 


1.640 


9.6108 


9.5379 


9.6581 




2S 


0.660 


0.560 


0.690 


9.7838 


9.7528 


9.8426 




IP 


0.578 


0.460 


0.640 


9.8226 


9.8274 


9.8670 



Table II shows that the critical values of the screening masses for cc and bh dissociations 
given by our calculations are larger than those given by Refs.[9,10]. This indicates that 
the results are dependent on the model. So a finer calculation of the screening masses for 
cc and hh dissociation is needed . Because the J /ijj suppression is related to the colour 
screening, the dissociation of J/^ is more interesting at finite temperatures. According 
to our calculations, the critical value of the screening mass for J/i/^ dissociation is about 
/Uc=0.900 GeV (the corresponding screening length is 0.219 fm). This information is 
probably useful to the study and observation of J/^ production in high energy collisions. 

In view of Table II, we would like to further note that the masses of all bound states 
are affected slightly by the change of ^. As shown by Refs.[9,10], with the increment of fi 
the masses of the cc bound states and those of the higher hh bound states decrease, while 
the T mass increases, which is different from ours. In our case, the masses of the higher cc 
and hh bound states decrease with /x, while the J /ip and T masses increase with ^u. This 
means that the positive string tension part of the potential dominates and is reduced as fi 
increases; only for the T and J /ip does the second term in the right hand side of Eq.(14) 
give the main contribution in the relativistic formalism. 

One can also find in Table II that the relativistic correction for cc bound states is 
larger than that of hh ones. This is not surprising. Even before the detailed numerical 
calculations, this qualitative conclusion can be reached based on the following reasonable 
physical consideration. 

As mentioned above, the Schrodinger equation is the non-relativistic limit of the BS 
equation. Usually, the b quark is heavy enough compared to Aqcd, and one can expect 
that the relativistic correction is small for hh bound states. Nevertheless, the mass of c 
quark , rric, is not much larger than Aqcd- The relativistic correction may be large in 
the case of cc bound states, which is corroborated by the detailed numerical calculations 
hsted in Table II. 

Finally, we check the sensitivity of our results with respect to the vector-scalar mixing 
parameter x appearing in the potential, the numerical results are listed in Table III. 

On the basis of the analysis of the numerical results (see Table III), we come to the 
following conclusions: 

For the systems containing two heavy quarks the sensitivity of the mass spectra to 
the Lorentz structure of the confining potential is rather moderate, especially in the hh 
system. Therefore, our results based on pure scalar confining potential are insensitive to 
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the particular Lorentz structure of confinement. 

Table III. The dependence of cc and 66 masses (GeV) on the mixing parameter x. 





states 


x = 0.0 


x = 0.1 


x = 0.3 


z = 0.5 


x = 0.7 


a; = 0.9 


x = 1.0 


cc 


IS 


3.067 


3.085 


3.120 


3.155 


3.189 


3.223 


3.240 




2S 


3.663 


3.686 


3.730 


3.773 


3.814 


3.854 


3.874 




3S 


4.019 


4.053 


4.118 


4.180 


4.238 


4.293 


4.320 




IP 


3.526 


3.552 


3.602 


3.650 


3.698 


3.745 


3.768 


66 


IS 


9.436 


9.440 


9.448 


9.455 


9.462 


9.470 


9.473 




2S 


10.013 


10.017 


10.027 


10.036 


10.047 


10.053 


10.058 




3S 


10.343 


10.350 


10.363 


10.376 


10.390 


10.402 


10.408 




IP 


9.901 


9.905 


9.915 


9.924 


9.933 


9.943 


9.947 



4 Discussions and Conclusions 

In summary, we have studied the mass spectra, root-mean-square radii and dissociation 
of cc and 66 in the relativistic formahsm. We would like to point out that a relativistic 
treatment of quark -antiquark bound states, by means of the reduced BS equation, does 
imply some improvement in the description of cc and 66 meson mass spectra. This indi- 
cates that the relativistic description is necessary for a finer calculation of J /ijj at finite 
temperatures. Because the nic is not heavy enough compared to Aqcd, the relativistic 
effect can not be neglected for the description of J/V' dissociation in hot matter. The 
critical values of the screening masses have also been calculated and compared with the 
previous results given by the non-relativistic models. The critical value of the screening 
mass for J/ip dissociation is ^c=0.600 GeV in Ref.[10] and /Uc=0.700 GeV in Ref.[9], re- 
spectively. In the present calculations, however, ;Uc=0.900 GeV, which is larger than those 
from non-relativistic models. This means that the magnitude of the screening masses is 
model dependent. In order to get a finer evaluation of the screening mass of J/ip, one 
must have a good knowledge the confining potential in addition to a relativistic descrip- 
tion. Therefore the study of quark confining potential is of importance and should be 
further proceeded. 

Stimulating discussions with Weiqing Zhao and Wanyun Zhao are gratefully acknowl- 
edged by one of the authors (H.Zong). We thank Baltin Reinhard for carefully reading 
our manuscript. 
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Figure Captions 

Fig.l The dissociation energies for cc bound states. The sohd hnes indicate our results, 
the dotted ones are from Ref.[10]. 



Fig. 2 The dissociation energies for hh bound states. The sohd hnes represent our results, 
the dotted ones are from Ref.[10]. 
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